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Abstract. We discuss the existence and uniqueness of discontinuous solutions 
to adjoint problems associated with nonlinear hyperbolic systems of conserva- 
tion laws. By generalizing the Haar method for Glimm-type approximations to 
hyperbolic systems, we establish that entropy solutions depend continuously 
upon their initial data in the natural norm. 



1. Introduction 

We discuss the existence and uniqueness of discontinuous solutions to adjoint 
problems associated with nonlinear hyperbolic systems of conservation laws. By 
generalizing the Haar method [30j to encompass discontinuous solutions generated 
by Glimm [26j and front-tracking |19L 1231 I54L [8] approximation methods, we 
establish that entropy solutions depend continuously upon their initial data in the 
natural norm. In order to implement the proposed linearization method, we 
analyze linear hyperbolic systems with discontinuous coefficients, possibly in a non- 
conservative form. Our analysis begins with the key observation that while entropy 
solutions, by definition, contain compressive shocks only, the averaged matrix 

'A{u,v)^ I Df{u + 0{v -u))de 
Jo 

(associated with the flux of the system / and two entropy solutions m, v) may con- 
tain compressive or undercompressive shocks, but no rarefaction-shocks. This is 
an essential observation since, otherwise, rarefaction shocks would be a source of 
non-uniqueness and instability. The proposed method rests on geometric proper- 
tics of A{u, v) and on the construction of a weighted norm which is determined 
during the evolution and takes into account wave cancellation effects along gen- 
eralized characteristics. This strategy was presented by the author in 1998 (in a 
lecture at the University of Wisconsin-Madison) and appeared in print in [33^ I39j . 
In independent work, another proof of this continuous dependence property for 
genuinely nonlinear systems, based on earlier work on scalar conservation laws by 
Liu and Yang f48l, was obtained simultaneously to |33L I39j by Bressan, Liu, and 
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Yang [14^ I49j . In the present review, we cover the recent results on non-genuinely 
nonhnear systems estabhshed in [40 1, 141] . 

2. Continuous dependence property via a linearization method 

We consider solutions with small total variation, constructed by Glimm or front 
tracking schemes, satisfying the strictly hyperbolic system 

(2.1) dtu + d^f{u) = 0, u^u{t,x) eR^. 

The flux f{u) e M.^ need not be genuinely nonlinear and we solely assume that, 
for all relevant values of u, the matrix Df(u) has distinct eigenvalues Xj{u) and 
basis of left- and right-eigenvectors Ij (u) and rj (u), respectively. Given two entropy 
solutions u,v we define 

ip := V — u 

and introduce the averaged matrix 

A:^A{u,v):=[ Df{du+{l-9)v)dd. 
Jo 

Clearly, if u,v satisfy (|2.1[) . then -0 is a solution to the following linear hyperbolic 
system 

(2.2) dt^P + d^iAi^) = 0. 

We observe that the stability property for the linear system, i.e. 

for a sufficiently large class of matrices A and solutions ■(/;, implies the continuous 
dependence property for the nonlinear system 

(2.3) ||(u-t;)(t)|Ui(K) < Wiu-vmh^R)- 

The main difficulties are, first, to identify a suitable class of linear systems and 
solutions and, second, to construct a weighted norm that decreases in time. Our 
objective in this short presentation is to briefly review our linearization method and 
present several new results; we refer to |40L 141] for further details. 

Throughout, it is convenient to work with piecewise constant data and solu- 
tions. We suppose that the given matrix- valued field A is strictly hyperbolic with 
distinct eigenvalues Xj and eigenvectors lj,rj. The points of continuity, jump, and 
interaction in A are denoted by 

IR+ X M = C{A)UJ{A)UJ{A), 

respectively. 

Definition 2.1. A j -discontinuity {t,x) G >J{A) propagating at the speed X is 
said to be 

- compressive if Xj > X > Xj , 

- slow undercompressive if X < min(A^, A^), 

- fast undercompressive if X > max(Aj, A^'"), or 

- a rarefaction-shock Xj < X < X~^ . 
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We use the following decomposition 



J{A) C{A) U S{A) U U 7^(A). 



We first point out that if A contains rarefaction-shocks then the uniqueness property 
and, consequently, the continuous dependence property fail, as illustrated by the 
following example. 

Example 2.2. Consider the case that the speed coefficient is a rarefaction shock 
in a scalar equation, i.e. 



Then, the corresponding Cauchy problem admits infinitely many weak solutions, 
since within the wedge \x/t\ < 1 the solution can not uniquely determined from the 
given initial data by the method of characteristics. 

Let us summarize key facts of the proposed linearization method: 

(1) Rarefaction shocks may lead to instability and non-uniqueness. 

(2) However, the entropy condition assumed by the weak solutions u, v to 
nonlinear hyperbolic systems implies that all shock waves in u, v are com- 
pressive and, more importantly, that the averaged matrix A{u, v) can not 
contain rarefaction shocks. 

(3) stability estimates for rarefaction-free systems are established by defin- 
ing functionals that are equivalent to the distance and are generated 
along the time-evolution by a constructive approach. 

(4) Furthermore, in the course of this analysis, pointwise convergence prop- 
erties of Glimm-type schemes are required to establish the stability of 
certain nonconservative products. 

Details can be found in [33l IM HE [HI HQI HI]- This method was first 
investigated, in the norm rather than the one, in the earlier papers [521 145) . 

Let us conclude this section by recalling the pointwise convergence property of 
Glimm-type methods, as established by Glimm and Lax |27| and DiPerna [231 124) . 
Sec also Dafermos |20L 121] , Liu [4 7) , and Bressan and LeFloch [13) . 

When applied to the nonlinear hyperbolic system ()2.ip with initial data of small 
total variation, Glimm solutions, say — u^{t,x), converge locally uniformly to 
the entropy solution u of the corresponding Cauchy problem, in the sense that the 
following two properties hold for all but countably many to and all xq: 

(1) If {tQ,X()) continuity point for u, then for all e > there exists a neigh- 
borhood J\f{tQ,XQ) of that point and a real Hq > such that for all h < 
and {t,x) e Af^{to,xo) 



0, ^ = 'ip{t,x) (^R, 



with 




X <Q, 
a; > 0. 



\u'\t,x) 



u{to,xo)\ + \u{t,x) - u{tQ,xo)\ < e. 



(2) If {tQ,XQ) a discontinuity point for u, then there exists a shock curve 
t y{t) defined near Iq with y(to) = and, for e > 0, there exist a 
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neighborhood J\f^ {xo,to), a real /ig > 0, and an approximate shock curve 
t y'^{t) such that for h < and [t,x) G N''^{to,xo) 

\y'^{t)-y{t)\<s, 

provided {t,y^{t)), (t,y{t)) e {to, xq), and for {t,x) G A/'^(io, a;o) 

\u''{t,x)-u{to,xo±)\<e, 

\uit,x) - uito,XQ±)\ < e, x^y(t). 

3. Linear and nonlinear scalar equations 

Linearization involving one solution. Consider the two scalar equations 

(3.1) dt^ + d^iUiu)^:)^0, dtu + d^fiu) = 0, 

where / : R ^ M is strictly convex. We emphasize that, when u — v, all shocks in 
a = a{u, u) = fu(u) are compressive. 

Following |37j , weak solutions tp are sought as measures in the spatial variables 
and p.ip is defined in the sense of Volpert |57j , which is equivalent to choosing the 
family of straightlines in Dal Maso-LeFloch-Murat's theory [22]. We impose initial 
data V'OjWq such that Uq € BV{M.), the space of functions with bounded variation, 
and V'o S -MtO^), the space of bounded measures. 

First, we consider the Riemann problem associated with data m;, Ur and ipijipr, 
respectively. We distinguish between two cases: either u is a shock (with speed 
denoted by A) or u is a rarefaction. 

Theorem 3.1 (Riemann problem [37\ ). (1) // the solution u is a shock, 

then the solution to the corresponding Riemann problem is given by 

^{t, x)^4'i + (A - ^i)H{x -\t)+t (C, - C/)<5,_it, 

with 

Cr (A - !u(uT))^r, C'l := (A - fu{ui))ipi. 

(2) // the solution u is a rarefaction, the Riemann problem admits infinitely 
many solutions, in particular (if^, g ]R) 

1p{t,x) (1 " H{x ~ tfu{ui)) + 1prH{x - tfu{Ur)) 

+ 5{x - tfu{Ur)) - (p*S{x - tfu{ui)). 

The following remarks are in order: 

(1) Case of a shock wave: 

- The Riemann problem admits a unique solution. 

- A Dirac mass in ip propagates along the shock trajectory and its 
strength grows linearly in time, and the solution is solely a bounded 
measure. 

- However, under the compatibility condition C/ = C,., the solution ip 
is more regular and is a function of bounded variation. 

(2) Case of a rarefaction wave: 

- The Riemann problem admits (at least) a one-parameter family of 
solutions. 

- The solutions tp contain two propagating Dirac masses and are bounded 
measures, only. 
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- Still, the i/)-equation does admit bounded variation solutions, ob- 
tained for instance by taking (/s* = in Theorem 13. II 
We now turn our attention to the Cauchy problem and we search for solution 
that are bounded measures in the space variable. 

Theorem 3.2 (Existence and uniqueness theory (37) ). Let u he an 
entropy solution with hounded variation to a conservation law with convex 
flux f. 

(1) For every initial data ipa £ i"'^(R) H L°°{R), the initial value prohlem for 
the equation 

dt^ + dMu) ^) = 
admits at least one solution ip £ L°°(K+, provided the product 

fu{u)'4> is understood as a Volpert product. 

(2) When ipo £ A4b(M.) and dxUo < Cq, the prohlem admits at most one such 
solution. 

The proof of the existence part uses Lax formula for the entropy solution u 
and allows us to write an explicit formula for ip by tracking forward the initial 
data ipo along the generalized characteristics associated with u. The uniqueness 
part is based on an L^-type contraction argument, which strongly uses the entropy 
condition satisfied at jump discontinuities. 

A large literature is available on linear hyperbolic equations. More general 
speed coefficients are now covered by Bouchut and James [7] and Popov [53] . The 
connection with the study of finite difference schemes was extensively investigated 
by Tadmor [56] . Important progress was also made on multidimensional equations 
by Ambrosio [l], Colombini, Crippa, and Ranch [18j : these latter papers impose 
a restriction on the divergence of the velocity field and do not cover the equations 
and solutions in Theorem 13.21 

Linearization involving two solutions. Dealing with two solutions is more 
delicate, since now the average speed does not contain only compressive shocks. 

Theorem 3.3 (Sharp stability property [41p . Consider the linear hyperholic 
equation 

atV + <9,(s^) = 0, 

where a :— a{u,v) and u,v are bounded variation entropy solutions of a scalar 
conservation law with general flux / : M ^ M. Then, for all BV solutions tp 
generated hy front tracking one has 

WmhHR) + D2(t) + D3(t) < ||^(0)|Ui(R) 

D2(i) / E |A^-5-||V'^|dr. 
•^0 era) 

D3(0:= E \a--M\a+-a.\\^p.\+ f \a-f'{u)\\i:\dV£. 

5(2)U.F(S) •'^ 

Here, V^f denotes the absolutely continuous part of the total variation measure 
of a. It should be observed that the term T)2{t) provides the strongest, quadratic 
decay and is associated with Lax shocks. The term Tisit) is cubic in nature and 
involves undercompressive waves S{a) U J- (a), only. 
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The above result immediately applies to = v — u and yields a sharp con- 
tinuous dependence property which generalizes the standard Kruzkov's contraction 
property. 

Earlier results about the continuous dependence property for scalar conserva- 
tion laws concerned convex Rux, only. 

(1) For approximate piecewise smooth solutions constructed by the Glimm 
scheme, Liu and Yang |48j introduced an explicit functional in order to 
control the norm between two solutions. 

(2) For exact solutions with bounded variation, Dafermos |21) first derived 
a sharp version using the method of generalized characteristics. Later, 
Goatin and LeFloch [28j covered a class of hyperbolic equations and de- 
veloped the technique based on compressive and undercompressive dis- 
continuities. 

In |41j . we show that, for scalar conservation laws with arbitrary flux- function, 
the averaged coefRcient a contains no rarefaction shocks. We construct a decreas- 
ing weighted norm for the solution -ip. We establish that the nonconservative 
products (especially the term in D3) arising in the sharp stability estimate are 
stable, thanks to the pointwise convergence property of the Glimm-type schemes 
mentioned earlier. Such argument were first developed by LeFloch and Liu [42j in 
their version of the Glimm scheme for nonconservative systems. (See Section |6l) 

For scalar conservation laws. Definition [6] simplifies into the following form. 

Definition 3.4. A shock wave at a point {t,x) e >J(a) associated with the 
speed a is said to be 

- compressive i/ a_ > a > a+ , 

- slow undercompressive ifa< min(a~, a"*"), 



- a rarefaction-shock if a < a < a"*" . 

The following classification provides us with a characterization of the nature 
of shocks in the averaged speed coefficient associated with two entropy solutions. 
We consider a scalar conservation law with flux / : M — > M and a propagating 
discontinuity (u-,u+) whose propagation speed is 



a(u_,M+):= / duf{su^ + {l — s)u+)ds. 
Jo 

Fix a constant u G R, representing the (constant) local value of another solution. 
Then, whether the discontinuity 



is a compressive, undercompressive, or rarefaction shock is uniquely determined by 
the sign of 






a_ :— a(ii_,u), x < a{u^^u^)t, 
CLj^ :—a(u^,v), x>a{u^,u^)t^ 



Q :— (u± — v) (m+ — w_) (a+ — a_). 

Specifically, we have: 

• If u± — V have the same sign, then 




<0, S{a), 
> 0, T(a). 
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If u± — V have opposite signs, then 

'>0, C{a), 



< 0, 7^(S) 



Next, by recaUing that sohitions under consideration also satisfy the entropy 
condition, we arrive at: 

Proposition 3.5 (Fundamental property for scalar equations). Let u, v be two 
entropy solutions (with arbitrarily large total variation) to a scalar conservation law 
with general flux. Then, an (entropy admissible) shock wave of m or v can not be 
a rarefaction-shock of the averaged speed a := a{u, v). 

4. Linear and nonlinear hyperbolic systems 

Linearization involving one solution. We begin with: 

Theorem 4.1 (Existence theory for the Cauchy problem |41) ). Given an en- 
tropy solution u (with small total variation) to a strictly hyperbolic system 

dtu + dxfiu) = 

and given an initial data i/jq G A1f,(R), the initial value problem associated with the 
linear hyperbolic system 

dt^b + d,(Dfiu)4')^0 
admits a (possibly non-unique) weak solution iplt) G A^b(M), satisfying on every 
compact time interval 

\\^{t)\\M,{R) < ||V'o||.M,(R)- 

For an earlier result in the genuinely nonlinear case, see Crasta and LeFloch 

The system under consideration above contains compressive shocks, only. Note 
also that solutions are understood in the sense of Volpert or, more generally, 
Dal Maso-LeFloch-Murat [22) . Our proof of Theorem 14.11 uses the property of 
pointwise convergence mentioned earlier in Section [2] and the following two lem- 
mas. 

Lemma 4.2 (Continuous part). There exists a time-dependent, bounded mea- 
sure fi = /i(i), supported on the (countable) set of jump points of x ^-^ u{t,x), such 
that for a.e. t>{) 

Dfiu'')^'' ^Dfiu)^ + tx 
in the weak-star sense of measures. 

Lemma 4.3 (Jump part). If t £ {t!\T ) ^ y^{t) is an approximate shock 
curve of a genuinely nonlinear family in , converging to some shock curve t G 
(T,T) ^y{t), then for a.e. te {T,T) 

lim^''({/(t)})=V.({y(t)}). 

A weaker result holds for non-genuinely nonlinear characteristic fields. We note 
that, for almost every time t and at every jump point x of the limiting function 
w, the left- and right-hand traces v± = v±{t,x) and u± — u±(t,x) of the limiting 
solutions w, V satisfy a relation 

w+ - w_ = 7 (u+ - u_) G R^, 
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where 7 — j{t, x) is a scalar. 

The Dirac mass solutions found in |37j are now refered to in the literature as 
"delta-shock waves" ; see Zhang and Zheng [58' , Danilov and Shelkovich [SF , and 
Mitrovic and Nedeljkov |51j . Furthermore, a large literature is available on solutions 
to linearized hyperbolic systems involving a single entropy solution, especially by 
Bressan and Colombo [10] , Bressan, Crasta, and Piccoli [HI [9] , Ancona and Marson 
[3], Bressan and Marson [15 , who investigated hyperbolic systems using the notion 
of "first-order tangent vector" associated with piecewise smooth solutions; these 
papers encompass general entropy solutions to Temple class systems, as well as 
genuinely nonlinear systems. 

Linearization involving two solutions. The linearization of systems in- 
volving two solutions is comparatively more challenging. As far as the continuous 
dependence property of entropy solutions is concerned, one main result achieved by 
the linearization method is as follows: 

Theorem 4.4 (Continuous dependence property [33 |,l39j ). Consider the Cauchy 
problem for a strictly hyperbolic system of conservation laws with genuinely nonlin- 
ear fields. Then, any two entropy solutions u, v with small total variation generated 
by front tracking depend continuously upon their initial data: 

||i;(t)-?/(i)||ii(R) < ||«(0)-u(0)|Ui(R). 

Other proofs of the continuous dependence property were simultaneously and 
independently proposed by Liu and Yang |49j (Glimm scheme and explicit func- 
tional), and by Bressan, Liu, and Yang [14| (front tracking and explicit functional). 
The latter proofs are reviewed in the textbooks 120^ I32j , while our proof via the 
linearization method included in the textbook |39j and was further developed in 
[291140]. 

Our proof of Theorem l4.4l is based a linearization approach and the construction 
of a weighted norm. A generalization to non-genuinely nonlinear systems of con- 
servation laws is presented in [41] and uses the general existence theory established 
in [Ml US, as well as in [5l[6], [50], [2l[4]. 

Let us sketch our proof of Theorem 14.41 We begin by introducing the charac- 
teristic components a = (q^) of a solution 

ip{t,x) =: ^aj{t,x)rj{t,x), {t,x) £ C{A) nC(V'), 

and the weighted norm 

Hit)\\w{t} / y2\a]{t,x)\wj(t,x)dx, 

where < Wmin < Wj{t,x) < Wmax- For simplicity, we restrict the presentation 
here to the case that A = A(u, v), the averaged matrix associated with two entropy 
solutions. 

It is not difficult to compute the time-derivative of the weighted norm. We sup- 
pose that the weight w formally solves the adjoint system. Then, for any piecewise 
constant -0 the weighted norm satisfies 

j^\\m\\Mt)<J2 E pj{t,x)wjit,x)+^+it,x)w+it,x) 

x€j{A{t)) 
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for all but finitely many t, with 

f3-{t,x) := (X{t,x)-Xj.{t,x)) \aj{t,x)\, 
P+it,x) := (\j+{t,x)-X{t,x)) \a+{t,x)\, 

which we call the characteristic flux of the solution tjj. 

Lemma 4.5 (Signs of the characteristic flux). Across each i-shock, one has for 

±/3f<0, j<i, 



and for j 



± /3f > 0, j> i, 

Pf < 0, A, Pf > 0, 



sgn{af) 



±Pt>0, Si, ±/3f<0, J'i. 

This leads us to the following conditions on the weights: 

• Case j ^ i: one always has one favorable sign and one unfavorable sign. 

• Case j = i: only rarefaction shocks correspond to two unfavorable signs ! 

Definition 4.6. Fix a small k > 0. Consider a solution ip = ip{t,x) together 
with its characteristic flux Pj. Then, at any i-shock the j -characteristic flux (1 < 
j < N) is said to be dominant if 

K\p-\>\?,+ -n.\\f37\ + \A+-A-\J2\(3-\. 

k 

The dominant components enjoy the following properties. 
Lemma 4.7 (Signs of the characteristic components). For all j ^ i 
sgn{a^) = sgn{a~), j dominant, 

while for j = i 

sgn{a^), CiUTZi and i dominant, 
—sgn{a~), SiUj^i and i dominant. 

Hence, we see that the change of sign of the characteristic component a is 
directly related to the nature (shock/rarefaction or undercompressive) of the dis- 
continuity in A. 

Proposition 4.8 (Fundamental property nonlinear hyperbolic systems). Given 
two entropy solutions u, v with small total variation to a strictly hyperbolic, gen- 
uinely nonlinear system and let A be the averaged matrix. Then, the matrix 
A(t,x) := A{u,v){t,x) may contain compressive and undercompressive shocks but 
can not contain rarefaction-shocks. 

The proof based on a monotonicity property for the eigenvalue Aj: 

For all u-,v varying in a small neighborhood of 0, the averaged speed 
Ai(-, v) is strictly monotone along the z-shock curve from u_. 
- Moreover, if the i-shock (u_, u+) satisfles Lax shock inequalities 

\i{u-) > \i{u-,u+) > Xi{u+), 
then the averaged speed satisfies 

Xi{u-,v) > Xi{u+,v). 
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Let US sketch the proof of these properties. The right-hand state u+ = w+(e) 
can be viewed as a function of u_ and a parameter e varying in the neighborhood 
of 0, with 

ii+(e) = u- +eri{u-) + O(e^). 

Then, we compute 

Xi{u+,v) - Xi{u^,v) = eViXi{u-,v) ■ n{u^) + O(e^) 

= e ViAi(u_,-u_) • r,{u-) + O(e^) + 0{e \v - -u_|) > 0. 

Using now that Ai(u,u) = Xi(u) and the symmetric property A(u^v) = A(v^u), we 
obtain 2 ViAi ■ = VAi • > 0. Provided jej + |u — u_| < (5i is sufficiently smah, 
we conclude that the function Ai(-, v) is strictly monotone along the shock curve. 

If the shock satisfies the entropy inequalities and the normalization 

VAi • Ti > is chosen, then e < and we conclude from the above calculation that 
the averaged speed decreases from -u_ to u+. 

5. Fluid dynamics equations 

Finally, we turn the discussion to the compressible fluid equations with general 
equations of state, in either the form of the Lagrangian p-system 

(5.1) dtui - dxU2 = 0, dtU2 + dxp{ui) = 0. 
with p'{ui) < 0, or equivalently in the form of the Euler equations 

(5.2) dtUi + dx{uiU2) ^ 0, dt{uiU2) + dx{ui{u2)^ + p{ui)) ^ 0, 
with nowp'(iii) > 0. 

Theorem 5.1 (Continuous dependence property for fluid dynamics [41] ) . Con- 
sider solutions generated by front tracking and with small total variation of the fluid 
dynamics equations (jS.ip and (j5.2p Then, any two entropy solutions u, v together 
with their approximations Uh,Vh satisfy for all t > 

\\vh{t) - Vh{t)\\mii) < \\vh{0) - Vh{0)\\LHm + o{h) 
and, in consequence, 

||z;(i)-w(t)|Ui(R) < ||«(0)-u(0)||ii(R). 
More precisely, we prove the following sharp continuous dependence estimate 
\\vit) - u(t)|Ui(R) + T>2it) + T>3{t) < \\viO) - M(0)||ii(R), 

where 



T>2{t) 




— U- \ dr, 



with a similar expression for D3(t). Similarly to the case of scalar equations, 
the terms D2(t) and Dsit) are associated with compressive and undercompressive 
shocks in both characteristic families, respectively. 

The main steps of the proof of Theorem 15.11 are as follows: 
(1) The existence part follows from Iguchi-LeFloch's theory of non-genuinely 
nonlinear systems |34| . 

- One first approximates the pressure function by a function with 
finitely many inflection points. 

- One constructs interaction functional adapted to the problem. 
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- One derives uniform estimates that are independent of the number 
of inflection points. 

(2) The continuous dependence part rehes on the following steps: 

- A classification of discontinuities based on the density/specific vol- 
ume variable is established. 

- We show that the averaged matrix A{u, v) contains no rarefaction 
shocks. 

- A weighted norm is obtained by a constructive method by solving 
an evolution equation. 

The key observation of the linearization method to apply to fluid dynamics 
equations is provided by: 

Proposition 5.2 (Fundamental property for fluid dynamics). Ifu,v are two 
entropy solutions to the fluid dynamics equations (jS.f |) or (15. 2p with general equa- 
tion of state for the pressure and for arbitrary large total variation, then the aver- 
aged matrix A(t,x) := A(u,v)(t,x) can not contain rarefaction- shocks. 



6. DLM theory of nonconservative hyperbolic systems 



To complete this presentation, we want to stress the importance of small-scale 
phenomena for formulating a well-posed hyperbolic theory. Recall the notion of 
family of paths in the sense of Dal Maso-LeFloch-Murat [22]. Let $ : [0, 1] 



pN 



pN 



be a map satisfying: 



(1) <&(•; u+) is a path connecting u_ to w+ 

$(0; M+) = u_, <f>(l; M+) = 

TV^[o,i]($(-;m_,u+)) < 

(2) <f> is Lipschitz continuous in the graph distance 

dist(^$(-; <&(•; ^ I"- - u'_\ + \u+ - u'_^\. 

DEFINITION 6.1. Given u S BV^(M,M^) and g a Borel function, there exists a 
unique measure, called a nonconservative product and denoted by 



giu) dxU 



that is uniquely defined by the two properties: 

(1) If B is a Borel subset ofC(u), then fi{B) :— g{u) dxU. 

(2) If x £ J{u) is a point of jump then, with u± := u±{x). 



[0,1] 



gm-:u^,u+))ds^i-;u-,u+). 



This definition provides us with a notion of weak solutions to, in particular, 
nonconservative systems, once a family of paths is prescribed. In the conservative 
case, our definition is independent of the paths and is consistent with the distribu- 
tional definition 

\7h{u)dxU = dxh(u). 
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Definition 6.2. A bounded variation function u is a weak solution in the DLM 
sense if for every test-function 9 



dtOu dxdt + 



A{u) dxU 



dt = 0. 



The Riemann problem for systems of conservation laws was first solved by Lax 
For nonconservative system we have the following generalization of Lax's 
theorem. 

Theorem 6.3 (Riemann problem [22| ). Given a nonconservative, strictly hy- 
perbolic, genuinely nonlinear system and a family of paths <&, the Riemann problem 
admits an entropy solution (in the DLM sense) satisfying Lax shock inequalities. 

Generalized Hugoniot jump relations for nonconservative systems read: 



^0 

Note that wave curves are only Lipschitz continuous at the origin, which is in 
contrast with Lax's standard regularity result. The regularity of the wave 
curves is also investigated in Hayes and LeFloch [31] , Bianchini and Bressan [5l [6] , 
Iguchi and LeFloch 34 , and Liu and Yang [50] . 

The initial value problem for nonconservative systems was solved by LeFloch 
and Liu. 

Theorem 6.4 (Existence theory for the Cauchy problem [42|). Consider so- 
lutions constructed by the Glimm scheme for nonconservative, strictly hyperbolic, 
genuinely nonlinear system. Then, w'' = u'^{t,x) have uniformly bounded total 
variation and converge to an entropy solution u in the DLM sense and for all but 
countably many times 

(6.1) \A{u^)d,^u^] {t) -^\A{u)d^v\ (t). 



Observe that the 'almost everywhere' convergence with respect to the Lebesgue 
measure is not strong enough to derive the nonlinear stability property (|6.ip and the 
pointwise convergence properties recalled at the end of Section[2]are necessary. It is 
interesting to observe similarities and differences between in carrying out the theory 
of weak solutions for conservative systems or for nonconservative systems. In their 
work on the vanishing viscosity method, Bianchini and Bressan [6] were able to cover 
both conservative and nonconservative systems by the same arguments. As far as 
Glimm-type methods are concerned, a major distinction between conservative and 
nonconservative systems must be pointed out, when non-GNL are allowed. The 
following superposition property, essential in the theory of conservative systems, 
does not extend to nonconservative systems: if three states ui,Um,Ur satisfy the 
two relations 

-A {Urn " Ul) + f{Um) - f {ui) = -A {Ur - Um) + /("r) - /("m) = 

for some speed A, then one also has 

-A {Ur - Ul) + f{Ur) - f{ui) = 0. 

For further material on nonconservative systems we refer to 
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